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ABSTRACT 

With the Akeno Giant Air Shower Array (AG ASA), 581 cosmic rays above 
lO^^eV, 47 above 4 x lO^^eV, and 7 above 10^°eV are observed until August 
1998. Arrival direction distribution of these extremely high energy cosmic rays 
has been studied. While no significant large-scale anisotropy is found on the 
celestial sphere, some interesting clusters of cosmic rays are observed. Above 
4 X lO^^eV, there are one triplet and three doublets within separation angle 
of 2.5° and the probability of observing these clusters by a chance coincidence 
under an isotropic distribution is smaller than 1 %. Especially the triplet is 
observed against expected 0.05 events. The cos(^gc) distribution expected from 
the Dark Matter Halo model fits the data as well as an isotropic distribution 
above 2 x lO^^eV and 4 x lO^^eV, but is a poorer fit than isotropy above 
lO^^eV. Arrival direction distribution of seven lO^^eV cosmic rays is consistent 
with that of lower energy cosmic rays and is uniform. Three of seven are 
members of doublets above about 4 x lO^^eV. 

Subject headings: cosmic rays — galaxies: general — large-scale structure of 
universe — Galaxy: halo 
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1. Introduction 

Investigation on anisotropy of extremely high energy cosmic rays is one of the most 
important aspects to reveal their origin. In energies ^ lO^^eV, cosmic rays slightly deflect 
in the galactic magnetic field if they are protons of galactic origin, so that one could observe 
the correlation of their arrival directions with the galactic structure. Especially in the 
highest observed energy range, correlation of cosmic rays with the local structure of galaxies 
may be expected if their origins are nearby astrophysical objects and the intergalactic 
magnetic field is less than 10~^ gauss. 



In the 1980's, Wdowczyk, Wolfendale and their collaborators ( [Wdowczyk and 



Wolfendale 1984] , [Szabelski, Wdowczyk and Wolfendale 19SB| ) have shown that excess of 
cosmic rays from the direction of the galactic plane increases systematically with energy 
until a little above lO^^eV, though the available data was not statistically enough at that 
time. Gillman and Watson (1993) have summarized anisotropics in right ascension and 
galactic latitude combining the Haverah Park data set with the data sets from the arrays at 



Volcano Ranch ([Linsley 1980| ), Sydney ( [Winn et al. 1986| ) and Yakutsk ([Efimov et al. 1986| ) 



No convincing anisotropics were observed; but large amplitude of the second harmonics 
at (4-8) X IQ-'^^eV was reported. Ivanov (1998) showed, with the Yakutsk data set, a 
north-south asymmetry in the galactic latitude distribution which is the southern excess 
with 3.5 a deviation from an isotropic distribution in (5 - 20) x lO^^eV. 

Recently, we have shown a significant anisotropy with first harmonic amplitude of ~ 
4 % in (0.8 - 2.0) x 10^'^eV, which corresponds to the chance probability of 0.2 % due 
to fluctuation of an isotropic distribution ( [Hayashida et al. 1998| ). This anisotropy shows 
broad cosmic-ray flow from the directions of the galactic center and the Cygnus regions. In 
the higher energies, no signiflcant large-scale anisotropy was found. Bird et al. (1998) have 
shown the galactic plane enhancement in the similar energy range. These experiments show 
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that significant fraction of cosmic rays around lO^^eV come from galactic sources. 

In the much higher energy range > 4 x lO^^eV, Stanev et al. (1995) have claimed 
that cosmic rays exhibit a correlation with the direction of the supergalactic plane and the 
magnitude of the observed excess is 2.5 - 2.8 cr in terms of Gaussian probabilities. Their 
result was mainly based on the Haverah Park data set. In the same energy range, such 
large-scale correlation with the supergalactic plane was not observed in the data sets of 
the AGASA ( [Hayashida et al. 1996|) , SUGAR ( [Kewley, Clay and Dawson 1996|) and Fly's 



Eye ( pird et al. 1998|) experiments. However, AGASA observed three pairs of cosmic rays 
above 4 x lO^^eV within a limited solid angle of the experimental accuracy and the chance 
probability is 2.9 % if cosmic rays distribute uniformly in the AGASA field of view. Two 
out of three are located nearly on the supergalactic plane. If cosmic rays in each of these 
pairs come from the same source, the detailed study on energy, arrival time and direction 
distribution of these clusters may bring information on their source and the intergalactic 
magnetic field QSigl and Lemoine 1998| , [Medina Tanco 1998| ). 



In the observed energy spectrum, there are two distinctive energies: E ~ lO^^eV and 
4 X lO^^eV. The former is the energy where the spectral slope changes ([Lawrence, Reid 



and Watson"T99l[ , [Efimov et al. 1991[ [Bird et al. 199^ [^oshida et al. 1995|, P^akeda et al 



1998 ). This is interpreted as transition from galactic to extragalactic origin. The latter is 



the energy where the GZK effect ( preisen 1966[ [Zatsepin and Kuz'min 1966[) , which is a 



series of energy loss through interaction with the cosmic microwave background photons, 
becomes important on their propagation from sources. It is important to study whether the 
arrival direction distribution of cosmic rays changes at these energies. 

Recent result of the AGASA energy spectrum shows the extension beyond the expected 
GZK cutoff ( jlakeda et al. 1998[ ). Since the distance to sources of cosmic rays above 



the expected GZK cutoff is limited to 50 Mpc ( [Hill and Schramm 1985[ , [Berezinsky and 
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Grigor'eva 1988, Yoshida and Teshima 1993 ), their arrival directions may be correlated 
with luminous matter distribution if they are astrophysical source origin such as hot spots 
of radio galaxies ([Biermann and Strittmatter 1987| , P?akahara 1990| , [Rachen and Biermann 



1993] , [Ostrowski 1998|) , active galactic nuclei ( piandford 1976| , [Lovelace 1976| , |Rees et aT 



1982|) , accretion flow to a cluster of galaxies (|Kang, Rachen and Biermann 1997|) , relativistic 



shocks in gamma-ray bursts ( |Vietri 1995| , |Waxmann 1995| ), and so on. There is another 
possibility that most energetic cosmic rays are generated through decay of supermassive 
"X" particles related to topological defects ( Bhattacharjee and Sigl 1998| , reference therein). 
In this case, arrival directions of most energetic cosmic rays are not necessarily associated 
with luminous matters. If such particles are the part of Dark Matter and are concentrated 
in the galactic halo, anisotropy associated with our galactic halo is expected ([Kuzmin and 



Rubakov 1997] , Perezinsky, Kachelriess and Vilenkin 19971) . 



In this paper, we first examine large-scale anisotropy in terms of various coordinates 
using the data set of the Akeno Giant Air Shower Array (AG AS A) until August 1998, 
including the old data set of the Akeno 20 km^ array (A20) before 1990. Then we search 
for the small-scale anisotropy above lO^^eV with the AGASA data set. 



2. Experiment 

The Akeno Observatory is situated at 138° 30' E and 35° 47' N. AGASA consists of 111 
surface detectors deployed over an area of about 100 km^, and has been in operation since 
1990 ( IChiba et al. 1992| , phoka et al. 1997|) . A20 is a prototype detector system of AGASA, 
operated from 1984 to 1990 ( [leshima et al. 1986| ), and is a part of AGASA after 1990. 

Each surface detector consists of plastic scintillators of 2.2 m^ area. The detectors are 
placed with a separation of about 1 km. They are controlled and operated from a central 
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computer through optical fiber network. Relative time difference among the detectors are 
measured with 40 nsec accuracy; all clocks at detector sites are synchronized to the central 
clock and signal-propagation time in cables and electronic devices are regularly measured at 
start of each run (twice a day). The details of the AGASA instruments have been described 
in Chiba et al. (1992) and Ohoka et al. (1997). 

The accuracy on determination of shower parameters are evaluated through the 
analysis of a large number of artificial events. These artificial events are generated with 
taking account of air shower features and fiuctuation determined experimentally. Figure |l] 
shows the accuracy on arrival direction determination for cosmic-ray induced air showers 
as a function of energies. The vertical axis denotes the opening angle A6 between input 
(simulated) and output (analyzed) arrival directions. The opening angles including 68 % 
and 90 % of data are plotted. By analyzing artificial events with the same algorithm used 
above, the accuracy on energy determination is estimated to be ± 30 % above lO^^eV 
( lYoshida et al. 1995|) . 



Table |l] lists the number of selected events, N(E), with zenith angles smaller than 45° 
and with core locations inside the array area. Events below lO^^eV are used only a reference 
analysis in this paper. The difference of N(E > 3.2 x lO^^eV) / N(E > lO^^eV) between 
A20 and AGASA arises from the difference of detection efficiency of each system. Seven 
events are observed above 10^°eV, including one event after Takeda et al. (1998). 

3. Results 

Figure |^(a) shows arrival directions of cosmic rays with energies above lO^^eV on the 
equatorial coordinates. Dots, open circles, and open squares represent cosmic rays with 
energies of (1-4) x lO^^eV, (4 - 10) x lO^^eV, and > 10^°eV, respectively. The shaded 
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regions indicate the celestial regions excluded in this paper due to the zenith angle cut of < 
45°. The galactic and supergalactic planes are drawn by the dashed lines. "GC" designates 
the galactic center. Figure H(b) shows arrival directions of cosmic rays only above 4 x 
lO^^eV on the galactic coordinates. Details of the cosmic rays above 4 x lO^^eV are listed 
in Table ^ 

3.1. Analysis in the Equatorial Coordinates 

3.1.1. Harmonic Analysis 

In order to search for cosmic ray anisotropy, it is required to compare observed and 
expected event frequencies at each region. An expected frequency is easily estimated as far 
as the exposure in each direction can be obtained; the uniformity of observation time on 
solar time for several years, which results in the uniform observation in right ascension, 
is expected for a surface array detection system operating in stable like AGASA. The 
fluctuation of the observation time on the local sidereal time is (0.2 ± 0.1) % which is small 
enough compared with anisotropy in this energy range, so that the exposure (observation 
time X collection area) in right ascension is quite uniform. 

Figure ^ shows results of the first (left) and second (right) harmonics in right ascension. 
The amplitude (top), the phase (middle), and the chance probability (bottom) are shown 
in each energy bin. In the top panels of the harmonic amplitude, the shaded region is 
expected from statistical fluctuation of an isotropic distribution with the chance probability 
larger than 10 %. No significant anisotropy above this level is found above 3.2 x lO^'^eV. 
This is consistent with our previous paper ( piayashida et al. 1998|) , in which zenith angles 



up to 60° were used. 



3.1.2. Declination Distribution 

Figure ^ shows the dechnation distribution of events above lO^^eV (hght shaded 
histogram) and 10^'' eV (dark shaded histogram). A sohd curve is a third order polynomial 
function fitted to the light shaded histogram. This curve is consistent with the zenith 
angle dependence of the AGASA exposure and considered to be the expected distribution 
if cosmic rays distribute isotropically on the celestial sphere. Since the trigger efficiency 
is independent of energy above lO^^eV and zenith angle less than 45°, this distribution is 
applied to in higher energies. Excess with 2.5 a deviation is found m. 5 = [30°, 40°] and this 
will be discussed later. 

3.2. Analysis in the Galactic and Supergalactic Coordinates 

3.2.1. Galactic and Supergalactic Plane Enhancement 

If cosmic rays have origin associating with nearby astrophysical objects, we may expect 
cosmic-ray anisotropy correlated with the galactic or supergalactic plane. Figure ^ shows 
the latitude distribution on the galactic (left) and supergalactic (right) coordinates in three 
energy ranges of (1-2) x lO^^eV (top), (2-4) x lO^^eV (middle), and > 4 x lO^^eV 
(bottom). A solid line in each panel indicates the cosmic-ray intensity expected from an 
isotropic distribution. In order to examine any preference for arrival directions along the 
galactic and supergalactic planes, the plane enhancement parameter fE introduced by 
Wdowczyk and Wolfendale (1984) was used. The fE value characterizes the anisotropy 
expressed by: 

Iobs{h)/Iexp{h) = (1 - Je) + 1.402 JE exp(-62), (1) 

where h is galactic or supergalactic latitude in radians, lohsip) and lexpip) are observed and 
expected intensities at latitude h. A positive fE value suggests a galactic or supergalactic 
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plane enhancement, fE = indicates that arrival direction distribution is isotropic, and 
a negative /e shows depression around the plane. Figure ^ shows the dependence of /e 
on the primary energy for the galactic (left) and supergalactic (right) coordinates. Some 
excess can be seen around the supergalactic plane in the seventh energy bin (log(E[eV]) = 
[19.1, 19.2]), where f§'^ = 0.36 ±0.15. In other energies, the arrival direction distribution is 
consistent with an isotropic distribution. 

3.2.2. 9gc Distribution 

Figure |^ shows the cos(6'gc) distribution, where Ogc is the opening angle between 
the cosmic-ray arrival direction and the galactic center direction, with energies above 
lO^^eV (top), 2 X lO^^eV (middle), and 4 x lO^^eV (bottom). Histograms are the observed 
distribution and the solid curves are expected from an isotropic distribution. The observed 
distribution is consistent with the solid curve in all energy ranges. The dashed and dotted 
curves are expected from the Dark Matter Halo model ( [Berezinsky and IVlikhailov 1998|) 



and will be discussed in Section 4.2. 



3.3. Significance Map of Cosmic-Ray Excess/Deficit 

There is no statistically significant large-scale anisotropy in the above one-dimensional 
analyses. Here, we search for two-dimensional anisotropy with taking account of the angular 
resolution event by event. 

Figures ^ and |^ show the contour maps of the cosmic-ray excess or deficit with respect 
to an isotropic distribution above lO^^eV and 4 x lO^^eV, respectively. A bright region 
indicates that the observed cosmic-ray intensity is larger than the expected intensity and 
a dark region shows a deficit region. For each observed event, we calculate a point spread 
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function which is assumed to be a norinahzed Gaussian probabihty distribution with a 
standard deviation of the angular resolution A9 obtained from Figure |l]. The probability 
densities of all events are folded into cells of 1° x 1° in the equatorial coordinates. At each 
cell, we sum up densities within 4.0° radius for Figure ^ and 2.5° for Figure p. These radii 
are obtained from a/2 x A^, and they would make excess regions clearer. The reference 
distribution is obtained from an isotropic distribution. In these figures, small statistics of 
observed and expected events result in bright regions at the lower and higher declination 
and hence bright spots below 5 = 0° are not significant. Two distinctive bright regions are 
found in Figure H, which are broader than the angular resolution. They are referred to as 
broad clusters, such as the BCl (20^50"^, 32°) and BC2 ( iMO'", 35°). The member events 
within 4° radius of BCl are listed in Table |^. Four brighter regions in the middle declination 
are found in Figure ^. the CI - C4 clusters which are noted in the eighth column of Table 
||. The CI - C3 clusters follow the notation used in our previous analysis ( [Hayashida et al. 



19961) . The C2 cluster is observed in both energy ranges. 



In Figure ||, the contour map has eight steps in [— 3cr, +3a]; lower two steps below 
—1.5(7 are absent. The significance of deviation from an isotropic distribution are estimated 
to be 2.4 a at the C2 cluster, 2.7 a at the BCl cluster, and 2.8 a at the BC2 cluster. The 
arrival directions of cosmic rays around the BCl cluster are shown in Figure |TD|(a), and a 
radius of each circle corresponds to the logarithm of its energy. Shaded circles have energies 
above lO^^eV and open circles below lO^^eV. Figure |10|(b) shows the arrival time - energy 
relation, and open circles denote members of the BCl cluster. The members of the BCl 
cluster have energies between lO^^eV and 2.5 x lO^^eV and no excess of cosmic rays are 
observed below lO^^eV around this direction. Five members of the BCl cluster are observed 
around MJD 50,000. This cluster is in the direction of a famous supernova remnant — 
the Cygnus Loop which extends about 3° around (20''50'", 30° 34'). The BC2 cluster is 
the broader cluster without a clear boundary. The BCl and BC2 clusters contribute the 
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excess around S = 35° shown in Figure ^. The C2 and BC2 clusters are located near the 
supergalactic plane and lead the largest /J*^ value in Section p. 2.1 . 



For small statistics of observed events, Figure ^ reflects the arrival directions of 
individual events (open squares and open circles in Figure |^). The brightest peak is at the 
C2 cluster where three cosmic rays are observed against expected 0.05 events. It is possible 
that some of these clusters are observed by a chance coincidence. It should be noted, 
however, that two of these clusters — the doublet (CI) including the AGASA highest 
energy event and the triplet (C2) — lie near the supergalactic plane, as pointed in our 
previous analysis ( [Hayashida et al. 19961 ). The arrival directions (left) and arrival time - 
energy relation (right) for the CI (top) and C2 (bottom) clusters are shown in Figure |rT[ 
A radius of each circle in the left panels corresponds to the logarithm of its energy, and 
open circles in the right panels denote members of the CI and C2 clusters. Around the C2 
cluster, several lower energy cosmic rays are observed very close to the C2 cluster. 

3.4. Cluster Analysis 

The threshold energy of 4 x lO^^eV is one distinctive energy where the GZK effect 
becomes large as mentioned in Section |I|. It is, however, quite important to examine what 
kind of dependence on threshold energy is operating. 

To begin with, we estimate the chance probability of observing one triplet and three 
doublets from 47 cosmic rays above 4 x lO^^eV. A cluster of cosmic rays is defined as 
follows: 

1. Define the i-th event; 

2. Count the number of events within a circle of radius 2.5° centered on the arrival 
direction of the i-th event; 
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3. If this number of events exceeds a certain threshold value Nth-, the i-th event is 
counted as a cluster. 



This procedure was repeated for total 47 events and then the total number of clusters with 
Nth was determined. The chance probability Pch of observing this number of clusters under 
an isotropic distribution is obtained from the distribution of the number of clusters using 
10,000 simulated data sets. These simulated data sets were also analyzed by the same 
procedure described above. Out of 10,000 simulations, 32 trials had equal or more doublets 
{Nth = 2) than the observed data set, so that Pch = 0.32%. And Pch = 0.87% for triplets 
{Nth = 3). 

Then, the energy dependence for observing (a) doublets and (b) triplets are estimated 
and the results are shown in Figure |12[ When a new cluster is added above a threshold 
energy, a histogram changes discontinuously at that energy. At the maximum threshold 
energy where the triplet is detected, we find Pch = 0.16 % in Figure |12|(b). The narrow 
peaks of Pch ^ 0.1 % above 4 x lO^^eV in Figures |12|(a) result from the CI, C3 and C4 
doublets, and another doublet C5 is found just below 4 x lO^^eV. Here, these chance 
probabilities are estimated without taking the degree of freedom on the threshold energy 
into account. However, the chance probabilities are smaller than 1 % and don't vary 
abruptly with energies above 4 x lO^^eV. This means that the threshold energy of 4 x 
lO^^eV for doublet and triplet in Figure § may indicate any critical energy, and suggests 
that their sources are not very far being different from those below this energy. 



3.5. lO^OeV Events 

Seven events have been observed with energies above 10^'^eV, and their energies 
and coordinates are also listed in Table 0. Their declination are near 6 ~ 20° while an 
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isotropic distribution is shown by the sohd curve in Figure ^. To check whether these seven 
events distribute isotropically or not, we compare celestial distribution of seven lO^^eV 
events with that for events between lO^^eV and 10^°eV in ten different coordinates. The 
Kolmogorov-Smirnov (KS) test ( Press et al. 19881) was used for avoiding any binning effect. 



The results are summarized in Table ^. The smallest KS probability in Table | is 2.5 % 
for the declination distribution; but this probability becomes larger using data set above 
6.3 X lO^^eV. One interesting feature is that five 10^°eV cosmic rays come from south-west 
of the AGASA array, where the strength of the geomagnetic field component which is 
perpendicular to an air shower axis is larger than the other directions ( |Stanev and Vankov 
pWfl) . 

4. Discussion 

4.1. Comparison with other experiments 

Above 3 X lO^'^eV, no large-scale anisotropy has been found with the harmonic analysis 
and the f^ fit. Gillman and Watson (1993) summarized the f^ values using the data sets 
obtained mainly from the Haverah Park experiment. They obtained no significant deviation 
from f^ = 0. The result from the Fly's Eye experiment ( [Bird et al. 1998| ) is consistent with 



an isotropic distribution of cosmic rays with E ^ lO^^eV. The analysis with the Yakutsk 
data set ( |Ivanov et al. 19971) shows no significant galactic plane enhancement above lO^^eV. 



The results from all experiments are consistent with this work on no correlation of cosmic 
rays above lO^^eV with the galactic plane. This may implicate the extragalactic origin of 
cosmic rays above lO^^eV if they are mostly protons. 

The BCl, BC2 and CI - C5 clusters are found with energies > lO^^eV or ;:^ 4 x lO^^eV. 
The C2 and BC2 clusters lead the small preference along the supergalactic plane in the 
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energy range of log(E[eV]) = [19.1, 19.2]. With the data sets of Haverah Park, Yakutsk, and 
Valcano Ranch ( |Uchihori et al. 1996| ) and AGASA, another triplet is found at the position 
of the CI cluster within experimental error box on arrival direction determination. This 
triplet at the CI cluster position includes the AGASA highest energy event and a lO^'^eV 
Haverah Park event. It should be noted that these triplets at the CI and C2 positions are 
close to the supergalactic plane. 

4.2. Correlation with Galactic Halo 

Kuzmin and Rubakov (1997) and Berezinsky et al. (1997) have suggested a cosmic-ray 
source model associated with Dark Matter distribution in our galactic halo. In this model, 
most energetic cosmic rays are generated through decay of supermassive particles which are 
trapped in the galactic halo and thus distribute symmetrically around the galactic center. 
The arrival directions of most energetic cosmic rays, therefore, exhibit anisotropy at the 
Earth ( [Berezinsky 1998" ). From recent studies by Berezinsky and Mikhailov (1998) and 



Medina Tanco and Watson (1998), a significant anisotropy would be expected in the first 
harmonics of right ascension distribution, the amplitude of 40 % at phase about 250°, which 
is independent of the ISO and NFW models of dark matter distribution in the galactic halo. 
The ISO and NFW models are described in Kravtsov et al. (1997) and Navarro, Frenk 
and White (1996), respectively. This expected anisotropy is consistent with the results of 
the harmonic analysis above 4 x lO^^eV as shown in Figure |^. However, this amplitude is 
explained with statistical fluctuation of an isotropic distribution. 

As shown by the dashed and dotted curves in Figure |^, the ISO and NFW models of 
Dark Matter distribution in the galactic halo lead excess toward the galactic center. Table 
H shows the reduced-x^ values of the observed cos(6'g'c) distribution with the isotropic, ISO 
and NFW models. Although the distribution expected from the ISO and NFW models 
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are quite different from the observed distribution in energies above lO^^eV, the reduced-x^ 
values are close to one another above 2 x lO^^eV and 4 x lO^^eV. Above 2 x lO^^eV, all 
three models are acceptable and it is hard to distinguish one from another. 

4.3. Correlation with Nearby Galaxies 

In Section |3.4|, we calculated the chance probability of observing clusters under an 



isotropic distribution. If cosmic rays are astrophysical source origin, the non-uniform 
distribution of galaxies or luminous matters should be taken into account, as claimed by 
Medina Tanco (1998). He calculated trajectories of cosmic rays above 4 x lO^^eV in the 
intergalactic magnetic field under the assumption that flux of cosmic rays is proportional to 
the local density of galaxies. The expected distribution of cosmic-ray intensity is no more 
uniform and this may result in a strong anisotropy. This is different from the results in 
this paper so that our estimation of the chance probability of observing clusters under an 
isotropic distribution is experimentally reliable. However, his calculation shows important 
results: the C2 cluster is on top of a maximum of the arrival probability for sources located 
between 20 and 50 Mpc; and the CI cluster locates on a high arrival probability region 
for sources at more than 50 Mpc. This suggests the possibility that the members of these 
clusters are generated at different sources. One need accumulate further statistics to 
make arrival direction, time and energy relation to be clear ( [Medina Tanco 1998| , ^igl and 



Lemoine 1998| ) to distinguish whether the members of clusters come from a single source or 



unrelated sources. 



-16- 



4.4. Correlation with the Known Astrophysical Objects 



As mentioned in Section O, the BCl cluster is in the direction of the Cygnus Loop 
(NGC6992/95). From the Hillas confinement condition of (magnetic field x size) for cosmic 
ray acceleration ( [Hillas 1984|) , the magnetic field in the shock of the Cygnus Loop is too 



small to accelerate cosmic rays up to lO^^eV. And the observed energy distribution and 
bunch of arrival time of the cluster members don't favor the diffusive shock acceleration. 
Another possible candidate is PSR 2053+36 with the period of 0.2215 sec and the magnetic 
field of about 3 x 10^^ gauss ( [Manchester and 'I'aylor 1981| ). It may be plausible that such 
highly magnetized pulsar has accelerated cosmic rays up to lO^^eV within a short time 
( punn and Ostriker 1969] , poldreich and Julian 196"9|) . It is highly desired to search for any 



signals from this direction in other energy range around MJD 50,000. 

For the CI - C5 clusters and 10^°eV cosmic rays, coincidence with known astrophysical 
objects are searched for from three catalogs which are the second EGRET catalog 
( [Thompson et al. 1995[ , [Thompson et al. 1996[) , the CfA redshift catalog ( [Huchra et al. 



[1995[) , and the eighth extragalactic redshift catalog ( [Veron-Cetty and Veron 1998[) . The 



selection criteria are the following: (i) the separation angles within 4.0° from a member 
of each cluster, and 2.5° for the lO^^eV cosmic ray; (ii) the redshift within 0.02. In the 
CfA catalog, only QSOs/AGNs are selected. Candidate objects are listed in Table ^. Out 
of these objects, Mrk 40 (VV 141, Arp 151) is an interacting galaxy and may be most 
interesting. It should be noted that Al-Dargazelli et al. (1996) claimed that nearby colliding 
galaxies are most favored as the sources of clusters (regions of excess events) defined by 
them using the world data available before 1996. 
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5. Summary 

In conclusion, there is no statistically significant large-scale anisotropy related to the 
galactic nor supergalactic plane. The slight supergalactic plane enhancement is observed 
just above lO^^eV and arises mainly from the BC2 and C2 clusters. Above 4 x lO^^eV, one 
triplet and three doublets are found and the probability of observing these clusters by a 
chance coincidence is smaller than 1 %. Especially the triplet is observed against expected 
0.05 events. Out of these clusters, the C2 (AGASA triplet) and CI (doublet including 
the AGASA highest energy event or triplet together with the Haverah Park 10^'^eV event) 
clusters are most interesting; they are triplets found in the world data sets and are located 
near the supergalactic plane. One should wait for the further high-rate observation to 
distinguish whether the members of clusters come from a single source or different sources. 
The cos(^Gc) distribution expected from the Dark Matter Halo model fits the data as 
well as an isotropic distribution above 2 x lO^^eV and 4 x lO^^eV, but is a poorer fit 
than isotropy above lO^^eV. The arrival direction distribution of the 10^°eV cosmic rays is 
consistent with that of cosmic rays with lower energies and is uniform. It is noteworthy 
that three of seven lO^'^eV cosmic rays are members of doublets. The BCl cluster is in the 
direction of the Cygnus Loop or PSR 2053+36 region. It is desirable to examine any signals 
from this direction in other energy band around MJD 50,000. We hope other experiments 
in TeV - PeV regions to explore the CI - C5 clusters and 10^°eV cosmic ray directions. 

We are grateful to Akeno-mura, Nirasaki-shi, Sudama-cho, Nagasaka-cho, Ohizumi- 
mura, Tokyo Electric Power Co. and Nihon Telegram and Telephone Co. for their kind 
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suggestions on the analysis of the Dark Matter Halo hypothesis. M.Takeda acknowledges 
the receipt of JSPS Research Fellowships. The authors thank Paul Sommers for his valuable 
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Fig. 1. — Accuracy on arrival direction determination. Closed and open circles are the 
opening angles encompassed 68 % and 90 % data. 

Fig. 2. — Arrival directions of cosmic rays with energies above lO^^^'^eV on the (a) equatorial 
and (b) galactic coordinates. Dots, open circles, and open squares represent cosmic rays with 
energies of (1 - 4) x lO^^eV, (4 - 10) x lO^^eV, and > 10^°eV, respectively. The galactic and 
supergalactic planes are shown by the dotted curves. "GC" designates the galactic center. 

Fig. 3. — Results of the harmonic analysis. (Top to bottom, the amplitude, the phases and 
the chance probabilities of the first (left) and second (right) harmonics.) 

Fig. 4. — Declination distribution of the observed cosmic rays. (Light shaded histogram: > 
lO^^eV. Dark shaded histogram: > 10^*^ eV, the right-had vertical axis should be referred.) 

Fig. 5. — Galactic (left) and supergalactic (right) latitude distribution. (Top: (1-2) x 
lO^^eV. Middle: (2 - 4) x lO^^eV. Bottom: > 4 x lO^^eV.) 

Fig. 6. — Dependence of the plane enhancement factor on the energy. (Left: for the galactic 
coordinates. Right: for the supergalactic coordinates) 

Fig. 7.— cos(^Gc) distribution. (Top: > lO^^eV. Middle: > 2 x lO^^eV. Bottom: > 
4 X lO^^eV.) Here, Oqc is the opening angle between the cosmic-ray direction and the 
galactic center direction, with energies > lO^^eV (top), 2 x lO^^eV (middle), and 4 x lO^^eV 
(bottom). The solid, dashed and dotted curves indicate the distribution expected for the 
isotropic, ISO and NFW models, respectively. 

Fig. 8. — Significance map of cosmic-ray excess/deficit above lO^^eV. The dashed and 
dash-dotted curve indicate the galactic and supergalactic plane, respectively. 

Fig. 9. — Significance map of cosmic-ray excess/deficit above 4 x lO^^eV. The dashed and 
dash-dotted curve indicate the galactic and supergalactic plane, respectively. 
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Fig. 10. — BCl cluster, (a) Arrival directions of cosmic rays around the BCl cluster. Radius 
of each circle corresponds to log(E[eV]), and shaded and open circles have energies above 
lO^^eV and between 3 x lO^^eV and lO^^eV, respectively, (b) Arrival time - energy relation. 
Open circles denote members of the BCl cluster and dots are cosmic rays near the BCl 
cluster. After the vertical dotted line, A20 is combined into AGASA. 

Fig. 11. — Arrival directions and arrival time - energy relation for the CI and C2 clusters. 
Here, cosmic rays above lO^^eV are plotted. (See also Figure 10) 

Fig. 12. — Energy dependence of the chance probability of observing (a) doublets and (b) 
triplets. 
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Table 1. Number of events in the data set. 



Array > lO^^eV > 4 x lO^^eV > 10^"eV 

A20 59 7 

AGASA 522 40 7 

Total 581 47 7 
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Table 2. AGASA events 4 x lO^^eV. 



Date 


Timc(JST) 


Energy 






Coordinates^ 




Note 








a. 


b 


l^ 


\P 




84/12/12 


14:18:02 


6.81 X lOi'^ eV 


22'' 


21"' 


38.4° 


93.3° 


-15.7° 




84/12/17 


10:28:16 


9.79 


18'' 


29" 


35.3° 


63.5° 


19.4° 




86/01/05 


19:31:03 


5.47 


4/. 


38" 


30.1° 


170.4° 


-11.2° 


C4 


86/10/23 


14:25:15 


6.22 


14'' 


02" 


49.9° 


96.8° 


63.4° 




87/11/26 


17:49:20 


4.82 


21'' 


57" 


27.6° 


82.1° 


-21.1° 




89/03/14 


02:45:39 


5.27 


13" 


48" 


34.7° 


68.3° 


75.6° 




89/08/16 


08:32:01 


4.07 


5" 


51" 


58.5° 


154.5° 


15.6° 




90/11/25 


11:05:39 


4.51 


16'' 


17" 


-7.2° 


6.1° 


29.6° 




91/04/03 


00:32:40 


5.09 


15" 


47" 


41.0° 


65.7° 


51.5° 




91/04/20 


08:24:49 


4.35 


18" 


59" 


47.8° 


77.9° 


18.4° 


C3 


91/05/31 


13:07:04 


5.53 


3" 


37" 


69.5° 


136.6° 


11.2° 




91/11/29 


14:53:03 


9.10 


19" 


06" 


77.2° 


108.8° 


25.6° 




91/12/10 


18:59:10 


4.24 


0" 


12" 


78.6° 


121.0° 


15.9° 




92/01/07 


03:16:49 


4.51 


9" 


36" 


38.6° 


184.3° 


48.0° 




92/01/24 


12:26:17 


4.88 


17" 


52" 


47.9° 


74.8° 


29.4° 




92/02/01 


17:20:52 


5.53 


0" 


34" 


17.7° 


117.2° 


-45.0° 




92/03/30 


03:05:30 


4.47 


17" 


03" 


31.4° 


53.6° 


35.6° 




92/08/01 


13:00:47 


5.50 


11" 


29" 


57.1° 


143.2° 


56.6° 


C2 


92/09/13 


08:59:44 


9.25 


6" 


44" 


34.9° 


180.5° 


13.9° 




93/01/12 


02:41:13 


10.12 


8" 


17" 


16.8° 


206.7° 


26.4° 




93/01/21 


07:58:06 


4.46 


13" 


55" 


59.8° 


108.8° 


55.5° 




93/04/22 


09:39:56 


4.42 


1" 


56" 


29.0° 


139.8° 


-31.7° 




93/06/12 


06:14:27 


6.49 


1" 


16" 


50.0° 


127.0° 


-12.7° 




93/12/03 


21:32:47 


21.3 


1" 


15" 


21.1° 


130.5° 


-41.4° 


CI 


94/07/06 


20:34:54 


13.4 


18" 


45" 


48.3° 


77.6° 


20.9° 


C3 



The celestial coordinates are based on the J2000.0 coordinates. 
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Table 2. AG AS A events above 4 x lO^^eV (continued). 



Date 


Timc(JST) 


Energy 






Coordinates 




Note 








a 


5 


^G 


6^ 




94/07/28 


08:23:37 


4.08 X lOi'^ eV 


Ah 


56"' 


18.0° 


182.8° 


-15.5° 




95/01/26 


03:27:16 


7.76 


11'' 


14" 


57.6° 


145.5° 


55.1° 


C2 


95/03/29 


06:12:27 


4.27 


17'' 


37" 


-1.6° 


22.8° 


15.7° 




95/04/04 


23:15:09 


5.79 


12'' 


52" 


30.6° 


117.5° 


86.5° 




95/10/29 


00:32:16 


5.07 


1" 


14" 


20.0° 


130.2° 


-42.5° 


CI 


95/11/15 


04:27:45 


4.89 


4/. 


41" 


29.9° 


171.1° 


-10.8° 


C4 


96/01/11 


09:01:21 


14.4 


16'' 


06" 


23.0° 


38.9° 


45.8° 


C5 


96/01/19 


21:46:12 


4.80 


3" 


52" 


27.1° 


165.4° 


-20.4° 




96/05/13 


00:07:48 


4.78 


17'' 


56" 


74.1° 


105.1° 


29.8° 




96/10/06 


13:36:43 


5.68 


13'' 


18" 


52.9° 


113.8° 


63.7° 




96/10/22 


15:24:10 


10.5 


19'' 


54" 


18.7° 


56.8° 


-4.8° 




96/11/12 


16:58:42 


7.46 


21" 


37" 


8.1° 


62.7° 


-31.3° 




96/12/08 


12:08:39 


4.30 


16'' 


31" 


34.6° 


56.2° 


42.8° 




96/12/24 


07:36:36 


4.97 


14" 


17" 


37.7° 


68.5° 


69.1° 




97/03/03 


07:17:44 


4.39 


19" 


37" 


71.1° 


103.0° 


21.9° 




97/03/30 


07:58:21 


15.0 


19" 


38" 


-5.8° 


33.1° 


-13.1° 




97/04/28 


13:46:18 


4.20 


2" 


18" 


13.8° 


152.9° 


-43.9° 




97/11/20 


07:23:25 


7.21 


11" 


09" 


41.8° 


171.2° 


64.6° 




98/02/06 


00:12:26 


4.11 


9" 


47" 


23.7° 


207.2° 


48.6° 




98/03/30 


08:17:26 


6.93 


17" 


16" 


56.3° 


84.5° 


35.3° 




98/04/04 


20:07:03 


5.35 


11" 


13" 


56.0° 


147.5° 


56.2° 


C2 


98/06/12 


06:43:49 


12.0 


23" 


16" 


12.3° 


89.5° 


-44.3° 




97/04/10 


02:48:48 


3.89 


15" 


58" 


23.7° 


39.1° 


47.8° 


C5 
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Table 3. Members of the clustering events above about 10 eV. 



Name 


Date 


Energy 






Coordinates^ 














a 


b 


l^ 


\P 


BCl 


95/10/09 


1.47 X 


lO^^eV 


20^ 


50™ 


30.8° 


73.9° 


- 8.2° 




95/11/23 


1.68 




20^ 


54™ 


34.2° 


77.1° 


- 6.8° 




95/07/18 


1.31 




20^ 


42™ 


33.2° 


74.8° 


- 5.5° 




95/09/24 


1.33 




20^ 


41'" 


34.1° 


75.4° 


- 4.8° 




91/07/02 


1.10 




20^ 


55™ 


35.1° 


77.9° 


- 6.4° 




96/08/02 


2.29 




20^ 


55™ 


32.4° 


75.9° 


- 8.1° 




97/05/28 


1.06 




20^ 


50™ 


34.7° 


77.1° 


- 5.9° 




97/06/20 


1.11 




21^ 


02™ 


33.7° 


77.8° 


- 8.4° 



BC2 1^*40™ 35° 134° -27° 



^The celestial coordinates are based on the J2000.0 coordinates. 
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Table 4. Kolmogorov-Smirnov test for celestial coordinates. 



KS-Probability 








KS-Probability 


0.268 


Zenith Angle (9) 






0.867 


0.202 


Declination (6) 






0.025 


0.085 


Ecliptic Latitude 






0.449 


0.182 


Galactic Latitute 


(6^) 




0.540 


0.654 


Supergalactic Latitude 


(6^G) 


0.167 



Azimuth Angle (0) 
Right Ascension (a) 
Ecliptic Longitude 
Galactic Longitude {l'-^) 
Supergalactic Longitude {l^'^) 



Table 5. Reduced-y^ values of the cosf^r^r*) distribution with three models. 



GC) 



> IQi^cV 



> 2 X IQi^eV 



> 4 X IQi^eV 



isotropic distribution 

ISO model 

NEW model 



2.0 
11.8 
10.0 



1.7 
2.2 
1.9 



1.8 
1.7 

1.6 



Table 6. Astrophysical objects near the AGASA events. 



Event ID 



Astrophysical object 



CI 
C2 

970330 (1.5 X lO^OcV) 



Mrk 359 (z = 0.017) 

NGG 3642 (z = 0.005), Mrk 40 (z = 0.02), Mrk 171 (z = 0.01) 

H 1934-063 (z = 0.011) 



-30- 



8^ 



CD go 

< 

D) 

< 

O) 

.£ 4 

c 

Q 
Q. 

O 



0' 



T \ 1 1 1 1 \ 1 1 1 1 \ 1 \ 1 1 1 r 




90% 



68% 



J \ I I I I \ I I I I \ I \ I I I L 



18.0 18.5 19.0 19.5 20.0 20.5 

Log(Energy[eV]) 



Fig. 1.- 



-31- 








360 





f 


r... 






~~_^° 





(b) 


S^ 


o\ ^g-^ 


/ / 


y o 
\A 


V 


o 


i. o ■ — \ 


^v+30 





A \ n 
/\ \ 
\ \ 


o \ 


3 ^ 


yi g ^ 


p 


^ 


\ 




1° 


\ \ ° 


o 
o 


O \ 


□ 

/ 
/ 


^ y 


\^ 


T V\ 


-..._: i 




o ./Z^^ 


^-30 




Six 


^ 


-"""■'^^O ° 









Fiff. 2. 



-32- 



E 
< 

"5 -1 
o 



First Harmonics 



19.0 19.5 

Log(Energy[eV]) 



-3 I 

18.5 



Second Harmonics 




19.0 19.5 20.0 

Log(Energy[eV]) 



19.0 19.5 20.0 

Log(Energy[eV]) 



5 I 

18.5 



19.0 19.5 20.0 

Log(Energy[eV]) 




19.0 19.5 20.0 

Log(Energy[eV]) 



19.0 19.5 

Log(Energy[eV]) 



Fiff. 3. 






Number of Events 



D 

o 
o 

q] 

CD 
CQ 



Ul 

o 



en 
o 



o 
o 




00 
00 



Ul 



Number of Events (E >_10 ^%V) 



-34- 



$ 40 



19.0 19.3 

10 eV -10 eV 




■50 50 

Galactic Latitude[deg] 



S 40 




-50 50 

Supergalactic Latitude[deg] 




-50 50 

Galactic Latitude[deg] 



E 10 

E 



19.3 19.6 

10 eV ~ 10 eV 




-50 50 

Supergalactic Latitude[deg] 



E 10 

E 



; > K^nv 








1 


- 


- 


p 


X 


- 


,-— 






\ _ 
















1 


" 



-50 50 

Galactic Latitude[deg] 



E 10 

E 



> 19.6 
- 10 eV 




-50 50 

Supergalactic Latitude[deg] 



Fig. 5. 



-35- 



1.0 



c 
a> 
E 0.0 

(D 

o 

c 
to 



-1.0 



18.5 



19.0 19.5 

Log(Energy[eV]) 



20.0 



1.0 



CO 



c 

o 

E 0.0 

o 

o 

c 

CO 



-1.0 



18.5 



19.0 19.5 

Log(Energy[eV]) 



20.0 



Fig. 6.- 



-36- 



100 



^—> 
c 

CD 
> 

LU 



CD 



CD 



CO 

^—> 
c 

CD 
> 

LU 



CD 



50 



c 20 
CD 
> 
LU 



10 



1 1 1 r 



10 eV 





H \ \ h 



H h 



> 19.3 
^ 10 eV 



H ^ h 



H ^ \ h 



H h 



19.6 



_ 10 eV 



■ . ^ 



J I L 



± 



J I I L. 



-1.0 



-0.5 0.0 

cosOgc) 



0.5 



1.0 



Fiff. 7.- 



-37- 




RA 



Fig. 8.- 







Fig. 9.- 



-38- 






' ' ' 


1 1 1 1 1 


_ 




(b) 




- 








_ 










O 20.0 






-- 














>. 








D) 
















Q) 








C 








LU 
















O) 
















_] 








19.5 


- 




- 






o 


' 






o 
o 


- 






CP 

° . _ • . s 


- 


19.0 




" 





, •, , 1 : .' 



46000 



48000 



50000 



MJD 



Fig. 10.- 



-39- 





F) 


^.-A 












( 


^ oo 


o 


Q 

c 


f 

n 


o C 


o 




V,-' 


( 





RA 






C1 




20.5 






III 


1 1 1 1 1 1 






o 


O 20.0 




- 






o 

c 

LU 




o 


_i 






19.5 




. 






. 


19.0 






III 


1 1 1 1 1 1 



46000 



48000 50000 

MJD 





C2 








20.5 












II 




- 


O 20.0 








- 










0) 

c 

LU 

R' 




o 





o_ 


_i 










19.5 








■ o _ 


19.0 


o 




^ o 




III 


r;""7"7"T""r"7"7" 



48000 50000 

MJD 



Fig. 11. 



-40- 



O 
O 



2 - 



(a) Doublet 



11] 



(b) Triplet 





19.4 19.5 19.7 19.4 19.5 19.7 

Log(Threshold Energy[eV]) 



19.9 



Fig. 12.- 



